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We report the results of a study of the exclusive charmless semileptonic decays, B + — > r]^l + v 
and B° — > 7r~£ + i/, undertaken with approximately 464 million _BB pairs collected at the T(4S') 
resonance with the BABAR detector. The analysis uses events in which the signal B decays are 
reconstructed with a loose neutrino reconstruction technique. We obtain partial branching fractions 
for B + — > r]£ + i/ and B° — } ir~l + u decays in three and 12 bins of q 2 , respectively, from which 
we extract the /+(q 2 ) form-factor shapes and the total branching fractions B(B + — > rjH + v) = 
(0.36 ± 0.05 stat ± 0.04 SHSt ) x 10" 4 and B(B° -»■ -K~i + v) = (1.42 ± 0.05 stat ± 0.07 itfst ) x 10" 4 . We 
also measure B(B + -> r)'£ + i>) = (0.24 ± 0.08 stat ± 0.03 S!/St ) x 10" 4 . We obtain values for the 
magnitude of the CKM matrix element \V u t\ using three different QCD calculations. 

PACS numbers: 13. 20. He, 12.15.Hh, 12.38.Qk, 14.40.Nd 



I. INTRODUCTION 

A precise measurement of the CKM matrix [l| ele- 
ment \V u b\ will constrain the description of weak inter- 
actions and CP violation in the Standard Model. The 
rate for exclusive charmless semileptonic decays involv- 
ing a scalar meson is proportional to \V u bf+{q 2 )\ 2 , where 
the form factor f+(q 2 ) depends on q 2 , the square of 
the momentum transferred to the lepton-neutrino pair. 
Values of f+{q 2 ) are given by unquenched Lattice QCD 
(LQCD) calculations 0,11) reliable only at large q 2 (> 16 
GeV 2 ), and by Light Cone Sum Rules (LCSR) calcula- 
tions [j, @ , based on approximations only valid at small 
q 2 (< 16 GeV 2 ). The value of \V u b\ can thus be deter- 
mined by the measurement of partial branching fractions 
of charmless semileptonic B decays. Extraction of the 
f+(q 2 ) form- factor shapes from exclusive decays @ such 
as B° -> ir~£ + v Q and B + -> t}^'H + v may be used to 
test theoretical calculations Q . The values of the branch- 
ing fractions (BF) of the B + — » rj^(. + v decays will also 
improve our knowledge of the composition of charmless 
semileptonic decays and help constrain the size of the 
gluonic singlet contribution to the form factors for these 
decays [f|. 

In this paper, we present measurements of the par- 
tial BFs AB(B+ -> i]t+v,q 2 ) and AB(B° -> Tr-t+v,q 2 ) 
in three and 12 bins of q 2 , respectively, as well as the 
total BFs for all three decay modes. Values of the 
total BFs were previously reported in Refs. @, l9rJl2j|. 
We use the values of AB(q 2 ) for the B° -> it~£ + v 
mode with form- factor calculations 2-4] to obtain val- 
ues of \V u b\. Values of \V u b\ have previously been ex- 
tracted from B° — > ir~l + v measurements by CLEO @, 
BABAR fH Hoi HI an d Belle [ll|. A verv recent measure- 
ment by BABAR [14[ will be discussed in Section VII . 
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II. DATA SAMPLE AND SIMULATION 

We use a sample of 464 million BB pairs correspond- 
ing to an integrated luminosity of 422.6 fb _1 collected 
at the T(4)S) resonance with the BABAR detector 15] at 
the PEP-II asymmetric-energy e + e~ storage rings and a 
sample of 44 fb _1 collected approximately 40 MeV be- 
low the T{AS) resonance (denoted "off-resonance data"). 
Detailed Monte Carlo (MC) simulations are used to op- 
timize the signal selections, to estimate the signal ef- 
ficiencies, and to obtain the shapes of the signal and 
background distributions. MC samples are generated for 
T(45) — > BB events, e + e~ — > uu / dd/ ss / cc/ t + t~ (con- 
tinuum) events, and dedicated BB samples containing 
B° — > n^£ + i' and B + —> i]^£ + v signal decays. The sig- 
nal MC events are produced with the FLATQ2 genera- 
tor [161 ] and are reweighted to reproduce the f+(q 2 , a, cb) 
Becirevic-Kaidalov (BK) parametrization [17| . where the 
values of the shape and normalization parameters, a and 
cb, are taken from Ref. Q- The BABAR detector's ac- 
ceptance and response are simulated using the GEANT4 
package (IsT ]. 



III. EVENT RECONSTRUCTION AND 
CANDIDATE SELECTION 

We reconstruct the B° — > ir~ t + v and B + — > r](''£ + v 
decays. The 77 meson is reconstructed in the rj — > 77 and 
77 — > TT + ir~ir° decay channels (combined BF of 62%) while 
the 77' is reconstructed in the 77' — > hyk + t{~ channel, fol- 
lowed by the 77 — > 77 decay (product BF of 17.5%) [l8l ]. 
The 77' — > p°7 decay channel suffers from large back- 
grounds and we do not consider it. We carry out an 
untagged analysis with a loose neutrino reconstruction 
technique Q, thereby obtaining a large candidate sam- 
ple. 

Event reconstruction with the BABAR detector is de- 
scribed in detail elsewhere [lf|. Electrons (muons) are 
identified by their characteristic shower signatures in 
the electromagnetic calorimeter (muon detector), while 
charged hadrons are identified using the Cherenkov de- 
tector and dE/dx measurements in the drift chamber. 
The average electron (muon) reconstruction efficiency 
is 93% (70%), while its misidentification probability is 
< 0.2% (< 1.5%). The neutrino four-momentum, P v = 
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(\Pmiss\iIhniss), is inferred from the difference between 
the momentum of the colliding-beam particles Pbeams 
and the vector sum of the momenta of all the parti- 
cles detected in a single event ptot-, such that p m %ss — 
Pbeams - Ptot- To evaluate Etou tne energy sum of all 
the particles, we assume zero mass for all neutrals since 
photons are difficult to disentangle from neutral hadrons 
and we take the mass given by the particle identification 
selectors for the charged particles. In this analysis, we 
calculate the momentum transfer as q 2 = (Pb — Pmeson) 2 
instead of q 2 = {Pi + P v ) 2 , where P B , P meson and Pi are 
the four-momenta of the B meson, of the 7r, 77 or rf me- 
son, and of the lepton, respectively. With this choice, 
the value of q 2 is unaffected by any misreconstruction of 
the rest of the event. Here Pb has an effective value. To 
estimate this value, we first combine the lepton with a 
7T, r\ or rf meson to form the so-called Y pseudoparticle. 
The angle, 9 by, between the Y and B momenta in the 
T(4S I ) frame, can be determined by assuming B — > Yv. 
In this frame, the Y momentum, the B momentum and 
the angle Oby define a cone with the Y momentum as 
its axis and where the true B momentum lies somewhere 
on the surface of the cone. The B rest frame is thus 
known up to an azimuthal angle <fi defined with respect 
to the Y momentum. The value of q 2 is then computed 
as the average of four q 2 values corresponding to four 
possible angles, <j>, <j> + n/2, 4> + n, 4> + 3tt/2 rad, where 
the angle <fi is chosen randomly and where the four val- 
ues of q 2 are weighted by the factor sin 2 9b , Ob being 
the angle between the B direction and the beam direc- 
tion in the T(45) frame [l9j]. We note that, 8 by being a 
real angle, | cos 9 by\ < 1. We correct for the reconstruc- 
tion effects on the q 2 resolution (0.51 GeV 2 ) by applying 
an unre gula rized unfolding algorithm to the measured q 2 
spectra [201 ] . 

The candidate selections are optimized to maximize 
the ratio S/yJ(S + B) in the MC simulation, where S 
is the number of signal events and B is the total num- 
ber of background events. The continuum background is 
suppressed by requiring the ratio of second to zeroth Fox- 
Wolfram moments |2l| to be smaller than 0.5. This back- 
ground is further suppressed for B° — > n~£ + i> by selec- 
tions on the number of charged particle tracks and neu- 
tral calorimeter clusters 22] that reject radiative Bhabha 
and converted photon processes. We ensure that the mo- 
menta of the lepton and meson candidates are kinemat- 
ically compatible with a real signal decay by requiring 
that a geometrical vertex fit of the tracks of the two par- 
ticles gives a y 2 probability greater than 1% and that 
their angles in the laboratory frame be between 0.41 and 
2.46 rad with respect to the e~-beam direction, the ac- 
ceptance of the detector. To avoid J/tp —> decays, 
we reject B° —> ir~ fj, + v candidates if the Y mass corre- 
sponds to the J/tp mass. The electron (muon) tracks are 
required to have momenta greater than 0.5 (1.0) GeV in 
the laboratory frame to reduce the number of misiden- 
tified leptons and secondary decays such as D — > Xlv, 
J/tp , t and kaon decays. Furthermore, the momenta of 



the lepton and the meson are restricted to enhance signal 
over background. We require the following: for B° — > 
Ti-i+v decays, |$| > 2.2 GeV or > 1.3 GeV or 
\p}\ + \p*J > 2.8 GeV; for B+ -> rji+u decays, \p* e \ > 2.1 
GeV or > 1.3 GeV or \p' e \ + \p* \ > 2.8 GeV; and for 
B+ -> ri'i+v decays, |$| > 2.0 GeV or > 1.65 GeV 
or 0.69 x \pl \ + > 2.4 GeV (all asterisked variables 

are in the center-of-mass frame). For the B + — s- rf'^l+v 
decays, we restrict the reconstructed masses of the rf 
and r\ to lie in the intervals 0.92 < m,< < 0.98 GeV and 
0.51 < m„ < 0.57 GeV. For these decays, we also re- 
ject events with q 2 higher than 16 GeV 2 , since the signal 
is dominated by background in that range. Most back- 
grounds are reduced by g 2 -dependent selections on the 
angle (cos 8 'thrust) between the thrust axes of the Y and 
of the rest of the event, on the polar angle (9 m i SS ) asso- 
ciated with Pmiss , on the invariant missing mass squared 
( m miss = E miss ~ Wmiss] 2 ) divided by twice the miss- 
ing energy (E mlss = E beams - E tot ), and on the angle 
(cos^) between the direction of the W boson (£ and v 
combined) in the rest frame of the B meson and the di- 
rection of the lepton in the rest frame of the W boson. 
The q 2 selections are shown in Fig. [1] and their effects 
illustrated in Fig. [2] for B° -> tt~£ + v decays. In Fig. H 
a single vertical line indicates a fixed cut; a set of two 
vertical lines represent a g 2 -dependent cut. The position 
of the two lines correspond to the minimum and maxi- 
mum values of the cut, as shown in Fig. [I] The functions 
describing the q 2 dependence are given in Tables IVllVIII 
of the Appendix for the three decays under study. For 
B + — > rjl + v decays, additional background is rejected by 
requiring that | cos 8y \ < 0.95, where 8y is the helicity 
angle of the rj meson [16| . 

The kinematic variables AE — (Pb ■ Pbeams — s/2)/ y/s 
and m ES = y/ (s/2 + p B ■ Pbeams) 2 1 E 2 eams - p B are 
used in a two-dimensional extended maximum-likelihood 
fit [23J to separate signal from background. Here, y/s is 
the center-of-mass energy of the colliding particles and 
Pb = Pmeson + Pi + Pv , m the laboratory frame. We 
only retain candidates with \AE\ < 1.0 GeV and tobs > 
5.19 GeV, thereby removing a region with large back- 
grounds from the fit. On average, fewer than 1.14 can- 
didates is observed per event. For events with multiple 
candidates, only the candidate with the largest value of 
cos 81 is kept. The signal event reconstruction efficiency 
varies between 8.3% and 14.6% for B° — > ■n~l + i>, and 
1.4% and 2.6% for B + — > -qt + v decays (77 channel), de- 
pending on q 2 . It is 0.6% for both B + — > rji + v (7r + 7T~7r° 
channel) and B + — > r\' l + v decays. 

IV. BACKGROUNDS AND SIGNAL 
EXTRACTION 

Backgrounds can be broadly grouped into three main 
categories: decays arising from b — > uiv transitions 
(other than the signal), decays in other BB events (ex- 
cluding b uiv) and decays in continuum events. For 
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FIG. 1: Distributions of the selection values in the signal re- 
gion for the q 2 -dependent variables used in the analysis of 
B° — } tt~£ + v decays. The vertical axis represents the selec- 
tion value for a given q 2 value. We reject an event when its 
value is in the shaded region. 
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FIG. 2: (color online) Distributions in the signal region for the 
g 2 -dependent selections used in the analysis of B° — ► ty~£ + u 
decays. The arrows indicate the rejected regions. All the 
selections have been applied except for the one of interest. In 
each panel, the signal area is scaled to the area of the total 
background. 



the B° — > tt~£ + v mode only, in which there are many 
events, each of the first two categories is further split 
into a background category where the pion and the lep- 
ton come from the decay of the same B ( "same B" cat- 
egory), and a background category where the pion and 
the lepton come from the decay of different B mesons 
("both B" category). 

Given the sufficient number of events for the ir£is decay 
mode, the data samples can be subdivided in 12 bins of q 2 
for the signal and two bins for each of the five background 
categories. Two bins are used for each background cat- 
egory since the background q 2 spectra are not that well 
known and need to be adjusted in the fit when the num- 
ber of events is sufficiently large to permit it. The q 2 
ranges of the background binning for the B° —> ir~ £ + v 
decays are: [0,18,26.4] GeV 2 for the b — > ulv same B 
category, [0,22,26.4] GeV 2 for the b ->• ulv both B cat- 
egory, [0,10,26.4] GeV 2 for the other BB same B cate- 
gory, [0,14,26.4] GeV 2 for the other BB both B category 
and [0,22,26.4] GeV 2 for the continuum category. In each 
case, the q 2 ranges of the two bins are chosen to contain a 
similar number of events. All the signal and background 
events are fitted simultaneously for a total of 22 param- 
eters. For the r\^lv modes, a smaller number of events 
leads us to restrict the signal and each of the three back- 
ground categories to a single q 2 bin except for the signal 
in the r\lv mode when 77 — » 77, which is investigated in 
three bins of q 2 . 

We use the AE-ities histograms, obtained from the 
MC simulation as two-dimensional probability density 
functions (PDFs), in our fit to the data to extract the 



yields of the signal and backgrounds as a function of 
q 2 . As an initial estimate, the MC continuum back- 
ground yield and q 2 -dependent shape are first normal- 
ized to match the yield and q 2 -dependent shape of the 
off-resonance data control sample. This results in a large 
statistical uncertainty due to the small number of events 
in the off-peak data. To improve the statistical preci- 
sion, the continuum background, initially normalized to 
the off-peak data, is allowed to vary in the fit to the data 
for the tt£v and 77^(77) modes where we have a large 
number of events. The fit result is compatible with the 
off-peak prediction within at most 1 standard deviation. 
Because of an insufficient number of events, the b — > ulv 
background is fixed in the fit for the rf''iv modes, and 
the continuum contribution is also fixed for the 77^i/(37r) 
and v( Iv modes. Whenever a background is not varied 
in the fit, it is fixed to the MC prediction except for 
the continuum background which is fixed to its normal- 
ized yield and q 2 -dependent shape using the off-resonanc 
data. The background parameters which are free in the 
fit require an adjustment of less than 10% with respect 
to the MC predictions. For illustration purposes only, we 
show in Fig. [3]Ai? and uies fit projections in the signal- 
enhanced region for B° — > tt~ £ + v decays in two ranges of 
q 2 corresponding to the sum of eight bins below and four 
bins above q 2 = 16 GeV 2 , respectively. More detailed AE 
and ities fit projections in each q 2 bin are also shown in 
Figs. [5] and M of the Appendix for the B° -!• tt~£ + l> de- 
cays. The data and the fit results are in good agreement. 
Fit projections for B + — s- rj('H + v, only available below 
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FIG. 3: (color online) Projections of the data and fit results 
for the B° — > tt~£ + v decays, in the signal-enhanced region: 
(a,b) tubs with -0.16 < AE < 0.20 GeV; and (c,d) AE with 
niES > 5.268 GeV. The distributions (a,c) and (b,d) are pro- 
jections for q 2 < 16 GeV 2 and for q 2 > 16 GeV 2 , respectively. 



TABLE I: Fitted yields in the full q 2 range for the signal 
and each background category, total number of MC and data 
events, and values of \ 2 f° r the overall fit region. 



Decay mode 




■q£ + v 




Signal 


11778 ± 435 


888 ± 98 


141 ± 46 


b — > u£v 


27793 ± 929 


2201 (.fixed) 


204(fixed) 


Other BB 


80185 ± 963 


17429 ± 247 


2660 ± 82 


Continuum 


27790 ± 814 


3435 ± 195 


517{fixed) 


MC events 


147546 ± 467 


23953 ± 183 


3522 ± 68 


Data events 


147529 ± 384 


23952 ± 155 


3517 ± 59 


X 2 /ndf 


411/386 


56/52 


19/17 



q 2 = 16 GeV 2 , are shown in Fig. |4] Table U gives the 
total fitted yields in the full q 2 range for the signal and 
each background category as well as the \ 2 values and 
degrees of freedom for the overall fit region. The yield 
values in the B + — > rjl + v column are the result of the fit 
to the combined 77 and 3-7T modes. 



V. SYSTEMATIC UNCERTAINTIES 

Systematic uncertainties on the values of the par- 
tial branching fractions, AB(q 2 ), and their correlations 
among the q 2 bins have been investigated. These uncer- 
tainties are estimated from the variations of the resulting 
partial BF values (or total BF values for B + — > rf ' l + v de- 
cays) when the data are reanalyzed with different simu- 
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FIG. 4: (color online) Projections of the data and fit results 
for the B + — > ry- 1 ' ' £ + v decays, in the signal-enhanced region: 
(a,b) rnss with -0.16 < AE < 0.20 GeV; and (c,d) AE 
with mES > 5.268 GeV. The distributions (a,c) and (b,d) 
are projections for the B + — > r]l + v and B + — >• rj'l + v decays, 
respectively, both for q 2 < 16 GeV 2 . 



lation parameters and reweightings. For each parameter, 
we use the full MC dataset to generate new AE-m es dis- 
tributions ("MC event samples") by varying randomly 
only the parameter of interest over a complete (> 3a) 
Gaussian distribution whose standard deviation is given 
by the uncertainty on the specific parameter under in- 
vestigation. One hundred such samples are generated for 
each parameter. Uncertainties due to B counting and fi- 
nal state radiation are estimated by generating only one 
sample. Each MC sample is analyzed the same way as 
real data to determine values of AB(q 2 ) (or total BF val- 
ues for B + — > rj ' l + v decays). The contribution of the 
parameter to the systematic uncertainty is given by the 
RMS value of the distribution of these values over the 
one hundred samples. 

The systematic uncertainties due to the imperfect de- 
scription of the detector in the simulation are computed 
by using the uncertainties, determined from control sam- 
ples, on the tracking efficiency of all charged particle 
tracks, on the particle identification efficiencies of signal 
candidate tracks, on the calorimeter efficiencies (varied 
separately for photons and K®), on the energy deposited 
in the calorimeter by mesons as well as on their pro- 
duction spectrum. The reconstruction of these neutral 
particles affects the analysis through the neutrino recon- 
struction. The uncertainties due to the generator-level 
inputs to the simulation are given by the uncertainties 
in the BFs of the background processes b — > ulv and 
b —> civ, in the BFs of the secondary decays producing 
leptons [11], and in the BFs of the T(45) -> BB de- 
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TABLE II: Values of signal yields, AB{q 2 ) and their relative uncertainties (%) for B° -+ 7T B + -> T]£ + u and B + -> f/^+i/ 

decays. 



Decay mode 






TT~£ + l> 






r) £ + v 


/7 2 rancrp i (~^p\/ ^\ 


n 2 ^1 9 
y xii 


f7 2 ^1 6 

y xiu 


n 2 >1 6 
y /iu 


full rP 1 fan <xp 
i_u.ii y i mi 


n 2 ^16 
y xiu 


a 2 ^~\ 6 
y ^ ±u 


ViplH 




8499 1 




1 1 777 Pi 
lit i i .\j 


887 Q 




DF J 


u.oo 


i no 


u.oo 




u.ou 


n 94 


Statistical error 


3.9 


3.7 


7.6 


3.5 


12.5 


32.8 


Detector effects 


3.1 


3.5 


6.1 


4.0 


8.0 


8.8 


Continuum bkg 


0.9 


0.8 


1.0 


0.7 


0.3 


7.1 


B -»■ X u £^ bkg 


2.0 


1.7 


4.2 


2.0 


7.6 


6.7 


S -)■ XJz/ bkg 


0.6 


0.7 


1.8 


1.0 


1.2 


2.6 


Other effects 


2.3 


2.2 


3.2 


2.3 


3.4 


4.6 


Total uncertainty 


5.9 


5.9 


11.3 


6.3 


17.0 


35.8 



cays @. The B —> Xlv form- factor uncertainties, where 
X = (tt, p, uj, , D, D*), are given by recent calculations 
or measurements [18j. The uncertainties in the heavy 
quark parameters used in the simulation of nonresonant 
b — > vlv events are given in Ref. [H[ ■ We assign an un- 
certainty of 20% [25[ to the final state radiation (FSR) 
corrections calculated by PHOTOS 26]. Finally, the un- 
certainties due to the modeling of the continuum are es- 
tablished by using the uncertainty in its q 2 distribution 
shape and, when the continuum background is fixed, the 
uncertainty in the total yield, both given by comparisons 
with the off-resonance data control sample. 

The list of all the systematic uncertainties, as well as 
their values for the partial and total BFs, are given in Ta- 
bles En]] and EX] of the Appendix. The item "Signal MC 
stat error" in these tables includes the systematic uncer- 
tainty due to the unfolding procedure. The correlation 
matrices obtained in the measurement of the partial BFs 
are presented in Tables [Xj IXII and IXIII A condensed ver- 
sion of all the uncertainties is given in Table IH together 
with signal yields and partial BFs in selected q 2 ranges. 
The values given for the B + — > r]l + v decays are those 
obtained from the fits to the distributions of the r\ — > 77 
and 77 — > 7T + 7r~7r° channels combined. The larger relative 
uncertainties occurring in bin 12 of Table IvTjTI are due to 
poorly reconstructed events, and to the small raw yield 
in that bin. The former arises from the presence of a 
large number of low momentum pions and a large back- 
ground. This makes it difficult to select the right pion 
and results in a larger absolute uncertainty on the fitted 
yield. The small yield leads to a fairly large unfolding 
correction in this bin and thus to a considerably reduced 
unfolded yield. On the other hand, the unfolding process 
increases the absolute uncertainty only slightly. The re- 
duced yield together with the larger absolute uncertainty 
lead to the larger relative uncertainties reported in the 
table. 



VI. RESULTS 

The partial BFs are calculated for B° —> tt~ l + v and 
B + — s- r\l + v decays using the unfolded signal yields, the 
signal efficiencies given by the simulation and the BFs 
B(T(45) -> B°B°) = 0.484 ± 0.006 and B(T(45) -> 
B+B~) = 0.516 ± 0.006 ||. We obtain the total BFs 
B(B° -> it-e+v) = (1.42±0.05 stat ±0.07 S3/st ) x HT 4 , 
B(B+ -> r}i + u) = (0.36 ± 0.05 stQt ± 0.04 syst ) x 10" 4 and 
B(B+ -> r)'£+v) = (0.24±0.08 stat ±0.03 syst ) x HT 4 . 
The BF value for B + — > rj'l + v has a significance of 3.2ct 
when we take into account only the statistical uncer- 
tainty [27| ■ Taking into account the effect of the system- 
atic uncertainty which increases the total uncertainty by 
about 8% leads to a reduced significance of 3. Oct. The 
BF value, obtained from a fit to the combined 77 and 
3-7T channels of the B + — > r\l + v decays, is in good agree- 
ment with the weighted average of the total BFs obtained 
separately for the 77 and 37r channels. Consistent re- 
sults are obtained when dividing the final data set into 
chronologically-ordered subsets, electron only and muon 
only subsets, modifying the q 2 or the AE and tues bin- 
nings, and varying the event selection requirements. 

The experimental AB(q 2 ) distributions are displayed 
in Fig.[5]for B° — > tt~(. + v decays and in Fig.[6]for B + — » 
rjl + v decays, together with theoretical predictions. To al- 
low a direct comparison with the theoretical predictions, 
which do not include FSR effects, the experimental dis- 
tributions in these figures have been obtained with the 
efficiency given by the ratio of q 2 unfolded events gener- 
ated after all the cuts with a simulation which includes 
FSR to the total number of events generated before any 
cut and with no FSR effects i.e. with PHOTOS switched 
off. We obtain the f+{q 2 ) shape from a fit to these 
distributions. The x 2 function minimized in the fit to 
the f+(q 2 ) shape uses the Boyd-Grinstein-Lebed (BGL) 
parametrization [28| consisting of a two-parameter poly- 
nomial expansion. For the B° —> n~l + v decays, the fit 
gives ai/a = —0.63 ± 0.29 and 02/^0 = —6.9 ± 1.7, 
with P(x 2 ) = 92.1% as well as a value of |V u6 / + (0)| = 
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FIG. 5: (color online) Partial AB(q 2 ) spectrum in 12 bins of 
q 2 for B° — > iv~£ + v decays. The data points are placed in the 
middle of each bin whose width is defined in Table IVIIII The 
smaller error bars are statistical only while the larger ones 
also include systematic uncertainties. The solid green and 
black curves show the result of the fit to the data of the 
BK [TtJ and BGL [28| parametrizations, respectively. The 
data are also compared to unquenched LQCD calculations 
(HPQCD [|, FNAL [|) and a LCSR calculation [|. 
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FIG. 6: (color online) Partial AB(q 2 ) spectrum in 3 bins of 
q 2 for B + — > r\i^v decays. The data points are placed in 
the middle of each bin whose width is defined in Table IIXI 
The smaller error bars are statistical only while the larger 
ones also include systematic uncertainties. The data are also 
compared to a LCSR calculation [jj. 



(8.5 ± O.Sstat ± 0-^syst) x 10 -4 from the fit extrapolated 
to q 2 = 0. This value can be used to predict rates of 
other decays such as B — > tttt [2i|. For completeness, we 
also show the fit to the BK parametrization [l7| . which 
gives a BK = 0.52 ± 0.04, with P( X 2 ) = 28.6%. 

The q 2 distribution extracted from our data is com- 
pared in Fig. [5] to the shape of the form factors obtained 
from the three theoretical calculations listed in Tabic Hill 
the one based on Light Cone Sum Rules [1] for q 2 < 12 
GeV 2 , and the two based on unquenched LQCD [HQ for 
q 2 > 16 GeV 2 . We first normalize the form-factor predic- 
tions to the experimental data by requiring the integrals 
of both to be the same over the q 2 ranges of validity given 
in Table Mil for each theoretical prediction. Considering 
only experimental uncertainties, we then calculate the x 2 
probabilities relative to the binned data result for various 
theoretical predictions. These are given in Table Mil for 
the B° — > 7r~ t + v decays. All three calculations are com- 
patible with the data. As shown in Fig. [5J a LCSR calcu- 
lation [j| is compatible with the data for the B + — > rjl + v 
decays. It should be noted that the theoretical curves 
in Fig. [5] have been extrapolated over the full q 2 range 
based on a parametrization obtained over their q 2 ranges 
of validity. These extended ranges are only meant to il- 
lustrate a possible extension of the present theoretical 
calculations. 

We extract a value of |V^f,| from the B° —> ir~t + v 
AB(q 2 ) distributions using the relation: \V u b\ = 



v/A6/(t s oAC), where t b ° = 1.525 ± 0.009 ps is the 
B° lifetime and A( = r/|I4&| 2 is the normalized partial 
decay rate predicted by the form-factor calculations @- 
4]. The quantities AB and A£ are restricted to the q 2 
ranges of validity given in Table Mil The values of A£ 
are independent of experimental data. The values of |Vu&| 
given in Table ITTTI range from (3.1 — 3.7) x 10~ 3 . A value 
of | 1 could not be obtained from the B + — > r)l + v de- 
cays because the required theoretical input A£ is not yet 
available. 



VII. COMBINED BABAR RESULTS 

At first glance, there appears to be a large overlap be- 
tween the present analysis of the B° — > ir~£ + i/ data and 
that of another recent BABAR measurement [lij]. How- 
ever, there are significant differences between the two 
analyses. Considering the same fit region, we obtain 
147529 selected events (signal or background) compared 
to 42516 such events in Ref. [l4j]. This difference can eas- 
ily be explained by the fact that we use the full BABAR 
data set in the present analysis but this is not so in 
Ref. Furthermore, the use of the loose neutrino 

reconstruction technique in this work leads to a larger 
background. Only 140 events are found in common be- 
tween the two data sets, i.e. 0.3% overlap. The statis- 
tical uncertainties are thus expected to be uncorrelated 
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TABLE III: Values of \V u b\ derived from the form-factor calculations for the B — > ty £ + v decays. The three uncertainties on 
\V u b\ are statistical, systematic and theoretical, respectively. 



q 2 (GeV 2 ) AB (1CT 4 ) AC (ps' 1 ) \V ub \ (lO' 3 ) x 2 /ndf Prob( X 2 ) 

HPQCD [2] > 16 0.33 ±0.03 ± 0.03 2.02 ±0.55 3.28 ±0.13 ±0.15t f 7 5.0/4 28.8% 

FNAL [3] > 16 0.33 ± 0.03 ± 0.03 2.21±°;g 3.14 ± 0.12 ± 0.14±^;|| 6.4/4 17.4% 

LCSR [4] < 12 0.84 ± 0.03 ± 0.04 4001^ 3.70 ± 0.07 ± 0.08jlg;|| 6.2/6 39.9% 



TABLE IV: Values of quantities of interest and their averages obtained in the study of B° — > n~ l + v decays. The third 
uncertainty, given for the average values, is due to the form-factor calculation. It is not shown for the individual determination 
of \V u b\- The results for ai/ao, 02/00 and |V^i,/+(0)| in the column titled "Average" are actually from a fit to the combined 
data, as discussed in the text. 



Present work Reference [141 



Total BF 

\V u b\HPQCD X 10 3 
\Vub\FNAL X 10 3 
\V u b\LCSR x 10 3 

|Kb/+(o)| x 10 4 

BGL 01/ao 
BGL 112/00 



1.42 ±0.05 ±0.07 
3.28 ±0.13 ±0.15 
3.14 ±0.12 ±0.14 
3.70 ±0.07 ±0.08 

8.5 ±0.3 ±0.2 
-0.63 ±0.27 ±0.10 
-6.9 ±1.3 ±1.1 



1.41 ±0.05 ±0.07 
3.21 ±0.13 ±0.12 
3.07 ±0.11 ±0.11 
3.78 ±0.08 ±0.10 

10.8 ±0.5 ±0.3 
-0.82 ±0.23 ±0.17 
-1.1 ± 1.6 ± 0.9 



Average 
1.42 ± 0.04 ± 0.07 
3.23 ± 0.09 ±0.13lgj 
3.09 ± 0.08 ±0.12^j 
3.72 ± 0.05 ± 0.09+o j 
9.4 ±0.3 ±0.3 
-0.79 ±0.14 ±0.14 
-4.4 ±0.8 ±0.9 



between the two analyses. The event reconstruction and 
simulation are also somewhat different. For example, the 
values of q 2 are computed using different, although in 
principle equivalent, relations: here, q 2 = (Pb — P-k) 2 
versus q 2 = (P t + P v f in Ref. Nevertheless, al- 

most all of the systematic uncertainties are expected to 
be highly correlated. 

It is gratifying to note that, as shown in Table IIV1 
the total BF as well as the values of \V u b\ obtained in 
the two analyses are in good agreement with each other. 
The value of |V u6 | quoted under Ref. 0] in Table [TV] for 
the FNAL [3J theoretical prediction is obtained using the 
values of the partial BFs given in Ref. [l4| for q 2 > 16 
GeV 2 . Similar numbers of signal events (11778 ± 435 
here compared with 10604 ± 376 in Ref. [Tj] when the 
events from B + — > ifil + v decays are also included) lead 
to similar statistical uncertainties in the two analyses. 

It is possible to obtain a good approximation to the 
average of the present results and those of Ref. [14| ob- 
tained in the B° — » tt~£ + i> decays by taking the statis- 
tical uncertainties to be uncorrelated and the systematic 
uncertainties to be fully correlated. Such an averaging 
procedure yields the averages, and associated uncertain- 
ties, given in Table lTVl for the total branching fraction and 
the values of \V u b\. The additional value of \V ub \ obtained 
in Ref. [l4{ with a combined fit to data and theoretical 
points is not taken into account in this computation. 

The above averaging method is not appropriate for the 
fitted BGL coefficients (ai/ao and 02/00) and the value of 
|T4b/+(0)|, since, as shown in Table HVl the two measure- 
ments of these quantities are only marginally compatible. 
Instead, we perform a new fit of the BGL parametrization 



to the combined partial branching fraction results from 
the two analyses, the 12 values obtained in this analysis 
and the six values from Ref. [T3|. Here again, the sta- 
tistical covariance matrices are uncorrelated and the sys- 
tematic covariance matrices are fully correlated between 
the two data sets. The combined error matrix from the 
two analyses is used to perform the fit, with the result 
shown in Fig. [7] and a \ 2 probability P(x 2 ) — 14.2%. 
When only the statistical covariance matrix is used, the 
X 2 probability is reduced to 3.1%. We note that the dis- 
crepancy in the two analyses of the partial BFs at low 
values of q 2 does not lead to discrepancies in the resulting 
values of the total BF or \V u b\, as is evident in Table ITVl 
Finally, we do not attempt to average the partial branch- 
ing fractions due to the different q 2 binning used in the 
two analyses. 



VIII. SUMMARY 

In summary, we have measured the partial BFs of 
B + — > rj£ + u decays in three bins of q 2 and of B° — s- 
iv~i + v decays in 12 bins of q 2 . From these distributions, 
we extract the f+(q ) shapes which are found to be com- 
patible with all three theoretical predictions considered 
for the B° — » it~ i + v decays and with the LCSR calcula- 
tion for the B + — > rj£ + v decays. The BGL parametriza- 
tion fits our data well and allows us to obtain the value 
of |T4b/+(0)|. Our measured branching fractions of the 
three decays reported in this work lead to a significant 
improvement in our knowledge of the composition of the 
inclusive charmless semileptonic decay rate. Our value 
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Unfolded q 2 (GeV 2 ) 



FIG. 7: (color online) Partial AB(q 2 ) spectrum for B° -> 
ir~£ + v decays, in 12 bins of q 2 for the present work and six 
bins of q 2 for Ref. [lj. The smaller error bars are statisti- 
cal only while the larger ones also include systematic uncer- 
tainties. The solid black curve shows the result of the fit to 
the combined data for the two analyses using the BGL [2f| 
parametrization. 



of the total BF for B + — > rf ' l + v is an order of mag- 
nitude smaller than the most recent CLEO result Q- 
Our value of the total BF for B + — > rj£ + v is consis- 
tent with a previous untagged BABAR result [l2|. The 
value of the ratio B(B+ -> ■q'£+v)/B(B + -> rj£ + v) = 
0.67±0.24 sta 4±0.11 sys t allows an important gluonic sin- 
glet contribution to the rf form factor. The present value 
of the total BF for B° — > tt~£ + v is in good agreement 
with a previous untagged BABAR measurement pj as well 
as with a recent BABAR result [l4| • It has comparable pre- 
cision to the present world average Q. For B° — > tt~ £ + v 
decays, we obtain values of \V u b\ for three different QCD 
calculations. The results are in good agreement with 
those of Refs. @, [HJ- The three values are all acceptable 



according to the data. Two of these values 0, H| are con- 
sistent, within large theoretical uncertainties, with the 
value measured in inclusive semileptonic B decays: \V u i>\ 
= (4.27 ± 0.38) x 1(T 3 Q. We also provide the aver- 
age values of the total BF and of \V„h\ obtained in the 
present work and those of Ref. 14]. We also give the 
values of |V u f,/+(0)|, <Xi/ao and 02/00 obtained in a com- 
bined BGL fit to the two data sets. It may be noted that 
our results for the B° ir~ £ + v decays are generally in 
good agreement with those obtained recently by the Belle 
collaboration [30l |. 
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IX. APPENDIX 

In Tables fVl I VII and I VIII we give the functions describ- 
ing the q 2 dependence of the selections used to reduce the 
backgrounds in the three decays under study. 

The list of all the systematic uncertainties, as well as 
their values for the partial and total BFs, are given in Ta- 
bles Em] and E] for the B° ^jrj+v and B+ ->■ rf'H+v 
decays, respectively. In Table IVIII1 we have one column 



for each bin of q 2 and three columns for various ranges of 
q 2 as well as the last column for the global result. In row 
1, "Fitted yield" , we give the raw fitted yield as the num- 
ber of events. In row 2, "Yield statistical error", we give 
the statistical uncertainty in % for each fitted yield. In 
row 3, "Efficiency" , we give the efficiency in % attached 
to each yield. In row 4, "Eff. (without FSR)", we give 
the efficiency in %, modified to remove the FSR effect. 
In row 5, "Unfolded yield", we give the yields from row 
1 unfolded to give the true values of the yields in each 
bin, expressed as the number of events. In row 6, 11 AB" , 
we give the values of the partial BFs computed as usual 
using the true (unfolded) yields and the efficiencies with 
FSR. In row 7, " AB (without FSR)", we give the values 
of the partial BFs computed as usual using the true (un- 
folded) yields and the efficiencies modified to remove the 
FSR effect. In rows 8 - 39, we give the contributions in 
% to the relative systematic uncertainties for each value 
of AB as a function of q 2 . In row 40, "Signal MC statis- 
tical error" , we give the statistical uncertainty due to the 
number of MC signal events. In row 41, "Total system- 
atic error" , we give the total systematic uncertainty in 
% for each value of AB, obtained as the sum in quadra- 
ture of all the systematic uncertainties in each column. 
In row 42, "Total statistical error", we give the statisti- 
cal uncertainty in % for each value of AB obtained from 
propagating the statistical uncertainties on the raw fitted 
yields, following the unfolding process and taking into ac- 
count the efficiencies. In row 43, "Total error", we first 
give the total uncertainty in % for each value of AB, ob- 
tained as the sum in quadrature of the total systematic 
error and the total statistical error. We then give, in the 
last four columns, the total uncertainties in % for each 
range of q 2 , obtained as the sum in quadrature of the 
total errors for the appropriate number of q 2 bins. A 
similar description applies to Table ITX1 

In our analysis, we compute the covariance matrix for 
each source of uncertainty, and use these matrices to cal- 
culate the uncertainties on the total BFs. The correlation 
matrices for the total statistical and systematic uncer- 
tainties are given in Table IXl for the B + — > i^ + v yields 
and in Tables El and EQ for the B° -> tt~£ + v yields, re- 
spectively. Finally, detailed AE and rriES fit projections 
in each q 2 bin are also shown in Figs. [5] and |H respec- 
tively, for the B° —> ir~ i + v decays. 
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TABLE V: g 2 -dependent selections used in B° — > ty^£ + u decays. 



cos9 e < 0.85 for all values of g 2 

cosO e > -0.0000167*g 8 +0.000462*g 6 +0.000656*g 4 -0.0701*g 2 -0.48 

mmiss/ZEmiss > —0.5 GeV for all values of g 2 

m 2 niss /2Emiss < 0.00544 * g 4 - 0.127 * g 2 + 1.37 GeV 

cos Qthruat < 0.9 for all values of g 2 

cos Othmst < -0.00159 * g 4 + 0.0451 * g 2 + 0.59 

Omiss > -0.000122 * g 6 + 0.00483 * g 4 - 0.0446 * g 2 + 0.405 rad 

(g 2 is given in units of GeV 2 ) 



TABLE VI: g 2 -dependent selections used in B + — > rj£ + u decays. 



cos 9t < 0.9 for all values of g 2 

cose e > 0.00629 * g 4 - 0.119 * g 2 - 0.252 

m 2 miss /2E miss < 0.8 GeV, g 2 < 7.5 GeV 2 

m 2 miss /2E miss < -0.05 * g 2 + 1.175 GeV, 7.5 < g 2 <16.0 GeV 2 

cose thruat < 0.05 * g 2 + 0.6, g 2 < 5.0 GeV 2 

cos dthvust < 0.85, 5.0 < g 2 < 16.0 GeV 2 

cos6 mlss < 0.92, g 2 < 11.0 GeV 2 

cose miss < 0.88, 11.0 < g 2 < 16.0 GeV 2 

(g is given in units of GeV^ 



TABLE VII: g 2 -dependent selections used in B + — > r/'£ + u decays. 



mmiss/ZEmiss > —0.3 GeV for all values of g 2 

m 2 miss /2E miss < 0.35 * g 2 + 0.325 GeV, g 2 < 2.5 GeV 2 

m 2 mias /2E misa < 1.2 GeV, 2.5 < g 2 < 4.5 GeV 2 

m 2 ni3S /2E miss < -0.1 * g 2 + 1.65 GeV, g 2 > 4.5 GeV 2 

cosOthrust < 0.05 * g 2 + 0.575, g 2 < 6.5 GeV 2 

cos6i thr „ st < 0.9, 6.5 < g 2 < 12.5 GeV 2 

cos Qthrust < -0.05 * g 2 + 1.525, g 2 > 12.5 GeV 2 

mi ss > -0.1 * g 2 + 0.45 rad, g 2 < 2.5 GeV 2 

6miss > 0.2 rad, 2.5 < g 2 < 5.5 GeV 2 

Omiss > 0.05 * g 2 - 0.075 rad, g 2 > 5.5 GeV 2 

(g 2 is given in units of GeV 2 ) 
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TABLE VIII: B° ->■ ix'tv yields, efficiencies (%), AZ3 (1(T 7 ) and their relative uncertainties (%). The AB and efficiency 
values labeled "without FSR" are modified to remove FSR effects. This procedure has no sigificant impact on the AB values. 



q' 2 bins ( GcV 2 ) 


0-2 


2-4 


4-6 


6-8 


8-10 


10-12 


12-14 


14-16 


16-18 


18-20 


20-22 


22-26.4 


g 2 <12 


g 2 <16 


g 2 >16 


Total 


Fitted yield 


894.7 


987.8 


1177.1 


1181.3 


1178.6 


1122.1 


996.1 


884.5 


904.3 


847.5 


729.9 


873.9 


6541.6 


8422.1 


3355.4 


11777.6 


Yield statistical error 


12.8 


8.1 


6.0 


6.4 


6.7 


7.0 


8.2 


9.8 


10.3 


10.5 


14.0 


21.0 


3.2 


3.6 


7.9 


3.7 


Efficiency 


8.34 


9.10 


9.22 


9.09 


8.59 


8.46 


8.53 


8.50 


9.40 


10.52 


11.61 


14.59 










Eff. (without FSR) 


8.00 


8.97 


9.15 


9.18 


8.63 


8.53 


8.58 


8.61 


9.45 


10.66 


11.71 


14.70 










Unfolded yield 


919.9 


960.7 


1189.6 


1184.5 


1182.9 


1141.5 


1027.3 


929.2 


979.5 


979.9 


905.8 


376.7 


6579.1 


8535.7 


3241.9 


11777.6 


AB 


122.7 


117.6 


143.6 


145.0 


153.4 


150.2 


134.1 


121.7 


116.0 


103.7 


86.8 


28.7 


832.5 


1088.3 


335.3 


1423.5 


AB (without FSR) 


128.0 


119.2 


144.6 


143.7 


152.5 


149.0 


133.3 


120.1 


115.3 


102.3 


86.1 


28.5 


837.1 


1090.5 


332.3 


1422.8 



Tracking efficiency 


3.2 


1.9 


3.1 


2.2 


2.3 


3.9 


2.6 


4.0 


3.5 


1.3 


4.1 


9.4 


2.3 


2.5 


2.9 


2.6 


Photon efficiency 


6.0 


3.4 


2.6 


1.3 


2.2 


2.5 


3.1 


3.0 


5.0 


1.4 


5.1 


24.2 


1.9 


2.2 


4.6 


2.7 


K ( l efficiency 


0.9 


0.3 


0.6 


0.3 


0.5 


0.4 


0.5 


0.8 


0.6 


0.4 


1.7 


6.8 


0.3 


0.3 


1.0 


0.4 


production spectrum 


1.0 


0.6 


1.0 


0.6 


1.1 


1.0 


0.6 


2.7 


1.7 


1.0 


2.0 


8.3 


0.7 


0.9 


1.9 


1.1 


Kt energy 

Li OJ 


1.0 


0.6 


0.3 


0.2 


0.2 


0.3 


0.2 


0.4 


0.6 


0.7 


0.8 


7.1 


0.2 


0.3 


0.8 


0.3 


£ identification 


4.0 


1.0 


1.2 


1.3 


0.6 


0.6 


1.6 


1.0 


0.9 


1.6 


0.7 


4.9 


0.3 


0.5 


1.1 


0.6 


7r identification 


0.5 


0.2 


0.2 


0.2 


0.2 


0.2 


0.2 


0.3 


0.3 


0.3 


0.3 


5.5 


0.2 


0.2 


0.7 


0.3 


Brcmsstrahlung 


0.6 


0.3 


0.1 


0.2 


0.2 


0.2 


0.2 


0.3 


0.3 


0.3 


0.3 


5.3 


0.2 


0.2 


0.7 


0.3 


g 2 continuum shape 


7.9 


1.6 


0.9 


0.3 


0.7 


0.3 


1.1 


0.3 


0.7 


1.1 


1.2 


5.4 


0.9 


0.8 


1.0 


0.7 


B(B+ -t TT°e+u) 


0.5 


0.1 


0.1 


0.1 


0.1 


0.1 


0.2 


0.3 


0.3 


0.3 


0.4 


7.1 


0.2 


0.2 


0.8 


0.3 


B(B a -t p-£+u) 


0.5 


0.3 


0.1 


0.1 


0.2 


0.1 


0.2 


0.3 


0.3 


0.3 


0.5 


10.1 


0.2 


0.2 


0.9 


0.3 


B(B+ -> p°£+u) 


0.7 


0.2 


0.1 


0.1 


0.2 


0.2 


0.2 


0.3 


0.4 


0.3 


0.3 


7.3 


0.2 


0.2 


0.8 


0.3 


B(B+ -> ul+u) 


0.5 


0.1 


0.1 


0.1 


0.1 


0.1 


0.2 


0.3 


0.3 


0.4 


0.3 


8.2 


0.2 


0.2 


0.9 


0.3 


B(B+ -> r)£+v) 


0.5 


0.1 


0.1 


0.1 


0.1 


0.1 


0.2 


0.3 


0.3 


0.3 


0.2 


5.5 


0.2 


0.2 


0.7 


0.3 


B(B+ -> ■q'i+v) 


0.5 


0.1 


0.1 


0.1 


0.1 


0.1 


0.2 


0.3 


0.3 


0.3 


0.3 


5.4 


0.2 


0.2 


0.7 


0.3 


Nonrcsonant b — > u£v BF 


0.6 


0.3 


0.1 


0.1 


0.2 


0.2 


0.2 


0.4 


0.5 


0.3 


0.6 


7.8 


0.2 


0.2 


0.7 


0.3 


SF parameters 


0.9 


0.5 


0.9 


0.4 


0.3 


0.5 


0.6 


0.3 


0.5 


2.3 


4.1 


23.4 


0.6 


0.4 


2.3 


0.8 


B -y p£v FF 


2.4 


1.4 


2.1 


1.4 


1.9 


1.6 


0.8 


0.7 


2.4 


3.0 


1.1 


16.6 


1.8 


1.5 


1.8 


1.5 


B° -> tt-£+u FF 


0.5 


0.1 


0.1 


0.1 


0.1 


0.1 


0.2 


0.3 


0.3 


0.3 


0.3 


7.5 


0.2 


0.2 


0.9 


0.3 


Other scalar FF 


1.1 


0.3 


0.5 


0.5 


0.4 


0.3 


0.3 


0.7 


1.7 


2.1 


2.0 


8.7 


0.4 


0.3 


0.6 


0.3 


B -> ui£v FF 


0.6 


0.1 


0.1 


0.1 


0.1 


0.1 


0.2 


0.3 


0.5 


0.6 


0.6 


18.1 


0.2 


0.2 


1.8 


0.5 


B(B -»■ D£u) 


0.6 


0.6 


0.2 


0.1 


0.3 


0.2 


0.5 


0.3 


0.4 


0.4 


0.4 


5.6 


0.2 


0.3 


0.7 


0.4 


B(B -> D*tv) 


0.7 


0.3 


0.2 


0.3 


0.4 


0.3 


0.3 


0.5 


0.4 


0.3 


0.5 


5.7 


0.3 


0.3 


0.7 


0.4 


B(B -> D**£v) 


0.7 


0.3 


0.3 


0.5 


0.5 


0.3 


0.8 


0.6 


1.1 


0.7 


0.7 


5.8 


0.3 


0.3 


0.9 


0.4 


Nonrcsonant b — > c£v BF 


0.6 


0.2 


0.2 


0.1 


0.2 


0.2 


0.2 


0.5 


0.3 


0.3 


0.3 


5.6 


0.2 


0.2 


0.7 


0.3 


B ->■ D^i/ FF 


0.5 


0.2 


0.3 


0.1 


0.2 


0.1 


0.2 


0.4 


0.4 


0.3 


0.4 


5.7 


0.2 


0.2 


0.7 


0.3 


S -»• FF 


0.6 


0.2 


0.2 


0.4 


0.2 


1.0 


0.6 


1.9 


0.4 


0.7 


1.1 


6.9 


0.2 


0.3 


1.0 


0.5 


T(45) -»■ B°B° BF 


1.5 


1.7 


1.2 


1.3 


1.3 


1.2 


1.5 


1.2 


1.4 


1.4 


1.0 


5.8 


1.3 


1.4 


1.2 


1.3 


Secondary lepton 


4.3 


3.2 


2.1 


1.2 


1.7 


0.5 


1.2 


0.5 


0.5 


0.9 


3.7 


5.8 


1.0 


0.9 


1.2 


0.9 


Final state radiation 


0.3 


1.3 


0.8 


2.2 


0.3 


1.4 


1.2 


1.3 


1.4 


1.6 


0.8 


3.4 


1.0 


1.1 


1.5 


1.2 


B counting 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


Fit bias 


0.1 


0.3 


0.4 


0.1 


0.1 


0.4 


0.4 


0.7 


0.1 


1.0 


2.0 


30.8 


0.2 


0.0 


1.8 


0.4 


Signal MC stat error 


1.3 


1.6 


1.4 


1.6 


1.4 


1.5 


1.4 


1.4 


1.3 


1.4 


1.2 


2.5 


0.6 


0.4 


0.6 


0.3 


Total systematic error 


12.9 


6.4 


6.0 


4.9 


5.0 


5.9 


5.6 


7.0 


7.8 


6.3 


10.0 


61.6 


4.4 


4.6 


8.3 


5.2 


Total statistical error 


14.7 


11.9 


9.0 


9.3 


9.4 


9.4 


10.8 


12.5 


12.8 


12.8 


17.6 


56.7 


3.9 


3.7 


7.6 


3.5 


Total error 


19.6 


13.5 


10.8 


10.5 


10.7 


11.1 


12.1 


14.3 


15.0 


14.3 


20.3 


83.8 


5.9 


5.9 


11.3 


6.3 



TABLE IX: B+ -> r) {,) l + v yields, efficiencies(%), AB (1(T 7 ) and their relative uncertainties (%). 



Decay mode 


If K. 1/ 






I ft V 


(77) 






(3-7T and 77 combined) 


q 2 bins ( GeV 2 ) 


Total 


Total 


0-4 


4-8 


8-16 


Total 


0-4 


4-0 


Q 1 fi 
O-10 


rcrtai 


Fittod viplH 


141.0 


244.8 


279.9 


216.8 


146.7 


643.4 


303.9 


331.5 


252.5 
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TABLE X: Correlation matrix of the partial A£>(£> + — > rjl + v,q 2 ) statistical and systematic uncertainties. 





Statistical 


Systematic 


q 2 bins ( GeV 2 ) 


0-4 4-8 8-16 


0-4 4-8 8-16 


0-4 


1.00 -0.08 0.00 


1.00 0.36 0.05 


4-8 


-0.08 1.00 -0.06 


0.36 1.00 0.29 


8-16 


0.00 -0.06 1.00 


0.05 0.29 1.00 
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TABLE XI: Correlation matrix of the partial AB(B° — ¥ n~£ + v,q 2 ) statistical uncertainties. 



q 2 bins ( GcV 2 ) 0-2 2-4 4-6 6-8 8-10 10-12 12-14 14-16 16-18 18-20 20-22 22-26.4 



0-2 1.00 -0.16 0.17 0.02 -0.02 0.03 0.01 0.04 0.05 0.02 0.04 -0.00 

2-4 -0.16 1.00 -0.32 0.11 0.00 -0.00 -0.01 0.01 0.01 -0.00 0.00 -0.00 

4-6 0.17 -0.32 1.00 -0.30 0.15 0.02 0.06 0.06 0.07 0.00 0.01 0.01 

6-8 0.02 0.11 -0.30 1.00 -0.22 0.13 0.07 0.06 0.07 0.00 0.00 0.02 

8-10 -0.02 0.00 0.15 -0.22 1.00 -0.22 0.16 0.05 0.08 0.01 -0.00 0.02 

10-12 0.03 -0.00 0.02 0.13 -0.22 1.00 -0.15 0.10 0.07 -0.01 0.02 0.00 

12-14 0.01 -0.01 0.06 0.07 0.16 -0.15 1.00 -0.16 0.13 -0.01 0.05 -0.00 

14-16 0.04 0.01 0.06 0.06 0.05 0.10 -0.16 1.00 -0.01 0.01 -0.02 -0.02 

16-18 0.05 0.01 0.07 0.07 0.08 0.07 0.13 -0.01 1.00 -0.17 0.09 -0.08 

18-20 0.02 -0.00 0.00 0.00 0.01 -0.01 -0.01 0.01 -0.17 1.00 0.05 -0.05 

20-22 0.04 0.00 0.01 0.00 -0.00 0.02 0.05 -0.02 0.09 0.05 1.00 -0.35 

22-26.4 -0.00 -0.00 0.01 0.02 0.02 0.00 -0.00 -0.02 -0.08 -0.05 -0.35 1.00 



TABLE XII: Correlation matrix of the partial AB(B° — > ty l + v,q 2 ) systematic uncertainties. 



q 2 bins ( GeV 2 ) 0-2 2-4 4-6 6-8 8-10 10-12 12-14 14-16 16-18 18-20 20-22 22-26.4 



0-2 1.00 -0.45 0.37 0.30 0.59 0.47 0.54 0.38 0.39 0.03 0.44 0.34 

2-4 -0.45 1.00 -0.24 0.03 -0.27 -0.09 -0.17 -0.19 -0.37 0.36 -0.37 -0.02 

4-6 0.37 -0.24 1.00 0.78 0.83 0.76 0.67 0.68 0.52 0.31 0.71 0.42 

6-8 0.30 0.03 0.78 1.00 0.71 0.74 0.65 0.63 0.46 0.47 0.53 0.32 

8-10 0.59 -0.27 0.83 0.71 1.00 0.74 0.71 0.70 0.50 0.35 0.66 0.38 

10-12 0.47 -0.09 0.76 0.74 0.74 1.00 0.74 0.80 0.61 0.36 0.61 0.38 

12-14 0.54 -0.17 0.67 0.65 0.71 0.74 1.00 0.69 0.73 0.29 0.56 0.33 

14-16 0.38 -0.19 0.68 0.63 0.70 0.80 0.69 1.00 0.71 0.34 0.65 0.36 

16-18 0.39 -0.37 0.52 0.46 0.50 0.61 0.73 0.71 1.00 -0.03 0.62 0.22 

18-20 0.03 0.36 0.31 0.47 0.35 0.36 0.29 0.34 -0.03 1.00 -0.02 0.18 

20-22 0.44 -0.37 0.71 0.53 0.66 0.61 0.56 0.65 0.62 -0.02 1.00 0.52 

22-26.4 0.34 -0.02 0.42 0.32 0.38 0.38 0.33 0.36 0.22 0.18 0.52 1.00 



0.5 " 1 
AE (GeV) 



05 1 
AE (GeV) 



0.5 " 1 "-1 
AE (GeV) 



8<q 2 <10GeV 2 



10<q 2 <12GeV 2 



12<q 2 <14GeV 2 



-0.5 0.5 1 

AE (GeV) 

14<q 2 <16GeV 2 




0.5 1 
AE (GeV) 



AE (GeV) 



0.5 1 -1 
AE (GeV) 



16<q 2 <18GeV 2 



18<q 2 <20 GeV 2 



20 < q 2 < 22 GeV 2 



-0.5 0.5 1 

AE (GeV) 

22 < q 2 < 26.4 GeV 2 




0.5 1 
AE (GeV) 



0.5 1 
AE (GeV) 



0.5 1 
AE (GeV) 



0.5 1 
AE (GeV) 



FIG. 8: (color online) AE yield fit projections in the signal-enhanced region, with tues > 5.2675 GeV, obtained in 12 q 2 bins from the fit to the experimental data for 
B° — > tt~£ + u decays. The fit was done using the full A.E-m,ES fit region. 
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FIG. 9: (color online) ttles yield fit projections in the signal-enhanced region, with —0.16 < AiJ < 0.2 GeV, obtained in 12 q 2 bins from the fit to the experimental 
data for B° — » ir~ £ + u decays. The fit was done using the full AE-tues fit region. 



